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ABSTRACT 

Characteristics  of  turbulence  in  the  lower  and  middle  troposphere  over  Owens  Valley  have  been  examined 
using  aircraft  in  situ  measurements  obtained  from  the  Terrain- Induced  Rotor  Experiment.  The  two  events 
analyzed  in  this  study  are  characterized  by  a  deep  turbulent  layer  from  the  valley  floor  up  to  the  midtropo¬ 
sphere  associated  with  the  interaction  between  trapped  waves  and  an  elevated  shear  layer.  Kelvin-Helmholtz 
(KH)  instability  develops  above  the  mountaintop  level  and  often  along  the  wave  crests  where  the  Richardson 
number  is  reduced.  The  turbulence  induced  by  KH  instability  is  characterized  by  a  progressive  downscale 
energy  cascade,  a  well-defined  inertial  subrange  up  to  1  km,  and  large  eddies  with  vertical  to  horizontal  aspect 
ratios  less  than  unity.  The  turbulence  below  the  mountaintop  level  is  largely  shear  induced,  associated  with 
wave  steepening  and  breaking,  and  is  more  isotropic.  Evaluation  of  structure  functions  indicates  that  while  the 
turbulence  energy  cascade  is  predominately  downscale,  upscale  energy  transfer  exists  with  horizontal  scales 
from  a  few  hundred  meters  to  a  few  kilometers  because  of  the  transient  energy  dispersion  of  large  eddies 
generated  by  KH  instability  and  gravity  wave  steepening  or  breaking. 


1.  Introduction 

Turbulence  is  ubiquitous  in  high  Reynolds  fluids  such 
as  the  earth’s  atmosphere.  Atmospheric  turbulence  has 
been  the  subject  of  numerous  studies,  mostly  based  on 
field  observations  or  wind  tunnel  experiments.  However, 
compared  to  boundary  layer  turbulence,  detailed  docu¬ 
mentation  of  atmospheric  turbulence  in  mountainous 
areas  is  still  relatively  rare.  Complexity  over  mountains 
is  enhanced  in  part  because  turbulence  can  be  generated 
by  mountain  waves  at  levels  well  above  the  atmospheric 
boundary  layer.  It  follows  that  observing  turbulence  over 
mountainous  areas  often  requires  relatively  expensive 
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platforms  such  as  well-equipped  research  aircraft.  Lilly 
and  Lester  (1974)  examined  stratospheric  waves  and  tur¬ 
bulence  using  high-altitude  aircraft  in  situ  measurements 
over  the  Sangre  de  Cristo  Mountains  of  Southern  Col¬ 
orado  obtained  from  the  Colorado  Lee  Wave  Obser¬ 
vational  Program.  They  found  that  multiple  peaks  were 
present  in  the  power  spectra  of  velocities  and  potential 
temperature,  associated  with  mountain  waves  and  wave- 
induced  turbulence.  The  mean  slopes  of  these  spectra 
were  approximately  -3  for  scales  between  3  and  20  km, 
likely  related  to  breaking  gravity  waves.  The  absence 
of  a  -5/3  slope  in  the  small-scale  limit  was  attributed 
to  observational  inaccuracy.  Strong  turbulence  was  ob¬ 
served  in  the  midtroposphere  over  central  Colorado  on 
11  January  1972  associated  with  a  severe  downslope 
windstorm  (Lilly  1978).  The  strong  vertical  wind  shear 
induced  by  large-amplitude  waves  and  the  subsequent 


©  2010  American  Meteorological  Society 


Form  Approved 
OMB  No.  0704-0188 


Report  Documentation  Page 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

JUL  2010 

2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2010  to  00-00-2010 

4.  TITLE  AND  SUBTITLE 

Turbulence  Characteristics  in  an  Elevated  Shear  Layer  over  Owens 

Valley 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

UCAR, Monterey, CA 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS (ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

Characteristics  of  turbulence  in  the  lower  and  middle  troposphere  over  Owens  Valley  have  been  examined 
using  aircraft  in  situ  measurements  obtained  from  the  Terrain-Induced  Rotor  Experiment.  The  two  events 
analyzed  in  this  study  are  characterized  by  a  deep  turbulent  layer  from  the  valley  floor  up  to  the 
midtroposphere  associated  with  the  interaction  between  trapped  waves  and  an  elevated  shear  layer. 
Kelvin?Helmholtz  (KH)  instability  develops  above  the  mountaintop  level  and  often  along  the  wave  crests 
where  the  Richardson  number  is  reduced.  The  turbulence  induced  by  KH  instability  is  characterized  by  a 
progressive  downscale  energy  cascade,  a  well-defined  inertial  subrange  up  to  1  km,  and  large  eddies  with 
vertical  to  horizontal  aspect  ratios  less  than  unity.  The  turbulence  below  the  mountaintop  level  is  largely 
shear  induced,  associated  with  wave  steepening  and  breaking,  and  is  more  isotropic.  Evaluation  of 
structure  functions  indicates  that  while  the  turbulence  energy  cascade  is  predominately  downscale,  upscale 
energy  transfer  exists  with  horizontal  scales  from  a  few  hundred  meters  to  a  few  kilometers  because  of  the 
transient  energy  dispersion  of  large  eddies  generated  by  KH  instability  and  gravity  wave  steepening  or 
breaking. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

18.  NUMBER 
OF  PAGES 

19a.  NAME  OF 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

17 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


2356 


JOURNAL  OF  THE  ATMOSPHERIC  SCIENCES 


Volume  67 


development  of  shear  instability  were  believed  to  be 
responsible  for  the  observed  midtropospheric  turbulence. 
Again,  the  -5/3  power  law  was  not  observed  in  the  ve¬ 
locity  and  potential  temperature  spectra.  Stratified  tur¬ 
bulence  associated  with  a  bora  event  was  examined  by 
Mahrt  and  Gamage  (1987),  and  a  variety  of  structure 
functions  were  evaluated  using  the  aircraft  data  ob¬ 
tained  from  the  Alpine  Experiment.  They  found  that  the 
observed  turbulence  was  strongly  influenced  by  stratifi¬ 
cation.  In  1999,  a  midtropospheric  mountain  wave  break¬ 
ing  event  over  the  Central  Alps  was  documented  during 
the  Mesoscale  Alpine  Program  (Jiang  and  Doyle  2004). 
The  velocity  spectra  derived  from  two  available  legs  ex¬ 
hibited  a  -5/3  slope  from  the  smallest  resolvable  scale 
up  to  approximately  4  km  and  the  mean  slopes  were  —2 
or  -3  for  larger  scales.  A  relatively  large  positive  buoy¬ 
ancy  production  rate  was  found  along  one  of  the  legs, 
implying  wave-induced  local  overturning  of  isentropic 
surfaces.  More  recently,  Doyle  et  al.  (2005)  studied  a  large- 
amplitude  mountain  wave  event  over  Greenland  observed 
during  the  Fronts  and  Atlantic  Storm-Track  Experiment 
(FASTEX).  Multiple  energy-containing  peaks  were  found 
from  the  velocity  spectra.  These  peaks  were  attributed 
to  large  eddies  induced  by  breaking  mountain  waves. 
Clearly,  understanding  turbulence  over  complex  terrain 
remains  a  challenging  problem.  In  mountainous  areas, 
turbulence  could  be  generated  by  multiple  mechanisms, 
including  wave-induced  shear  and  convective  instability 
associated  with  overturning  of  isentropic  surfaces,  and, 
in  the  lower  troposphere,  terrain-scale  circulations  driven 
by  differential  heating. 

The  Terrain-Induced  Rotor  Experiment  (T-REX)  took 
place  in  March  and  April  2006  over  and  within  Owens 
Valley  in  California.  Owens  Valley  is  a  quasi-two- 
dimensional  valley  located  between  two  nearly  parallel 
mountains,  the  Sierra  Nevada  and  the  Inyo  Mountains 
(Fig.  1).  Strong  mountain  waves  and  rotors  can  be  gen¬ 
erated  over  Owens  Valley  as  the  westerly  or  southwest¬ 
erly  flow  impinges  on  the  southeast-northwest-oriented 
topography.  The  objectives  of  T-REX  include  advanc¬ 
ing  our  understanding  of  mountain  wave-induced  rotors, 
interaction  between  waves  and  turbulence,  and  differen¬ 
tial  heating-forced  circulations  within  the  valley  (Grubisic 
et  al.  2008).  In  addition  to  the  dense  ground-based  in¬ 
struments  deployed  in  Owens  Valley,  the  University  of 
Wyoming  King- Air  (UWKA)  flew  stacked  legs  in  the 
lower  and  middle  troposphere  along  a  track  oriented 
approximately  across  Owens  Valley  (Fig.  1)  to  sample 
gravity  waves  and  turbulence  during  each  intensive  ob¬ 
servational  period  (IOP).  Approximately  two  dozen 
UWKA  flights  were  conducted  during  the  two-month¬ 
long  field  observational  period,  which  provided  a  valuable 
dataset  for  studying  waves  and  turbulence  over  complex 


Fig.  1.  The  terrain  height  in  the  observational  area,  including 
Owens  Valley,  central  Sierra  Nevada,  and  Inyo  Mountains,  is  shown 
in  grayscale  (interval:  0.3  km)  and  contours  (interval:  0.6  km).  The 
NCAR ISS  site  is  marked  by  a  triangle  and  the  UWKA  flight  track 
is  indicated  by  a  bold  line. 


terrain.  In  this  study,  we  analyze  two  UWKA  flights  with 
a  focus  on  the  characteristics  of  turbulence  over  Owens 
Valley.  In  particular,  we  want  to  deepen  our  understand¬ 
ing  of  issues  such  as  wave  and  turbulence  interaction, 
characteristics  of  large  eddies  associated  with  wave- 
induced  Kelvin-Helmholtz  (KH)  instability,  the  tur¬ 
bulence  inertial  subrange,  turbulence  intermittency,  and 
cross-scale  energy  transfer  over  complex  terrain. 

The  remainder  of  this  paper  is  organized  as  follows. 
The  relevant  observational  details  and  flow  conditions 
are  described  in  section  2.  General  wave,  eddy,  and  tur¬ 
bulence  characteristics  are  illustrated  in  section  3.  Break¬ 
ing  KH  billows  observed  along  four  midtroposphere  legs 
are  examined  in  section  4.  In  section  5,  the  turbulence 
characteristics  are  analyzed.  The  in  situ  measurements 
from  another  IOP  are  examined  in  section  6.  The  results 
are  summarized  in  section  7. 


2.  Background 

a.  King- Air  measurements 

During  each  westerly  IOP  (i.e.,  each  IOP  with  pre¬ 
vailing  westerly  or  southwesterly  winds  at  the  moun- 
taintop  level),  the  UWKA  flew  stacked  legs  along  a  track 
oriented  across  Owens  Valley  and  approximately  parallel 
to  the  wind  direction  at  the  mountaintop  level  (~4  km 
MSL).  The  altitude  of  flight  legs  ranges  from  1400  m  (i.e., 
approximately  300  m  above  the  valley  floor)  to  8500  m 
MSL  (Table  1).  The  leg  lengths  vary  approximately  from 
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Table  1.  IOP  4a  flight  traverse  summary.  The  possible  ranges 
for  upscale  energy  transfer  are  derived  from  third-order  structure 
functions. 


Time 

Data 

Upscale 

Leg  (UTC) 

Alt  (m)  Char  points  a 

range  (m) 

1 

1656-1702 

7590 

Wave 

9025 

0.10 

— 

2 

1705-1721 

5750 

Wave 

24  025 

0.12 

— 

3 

1723-1730 

6660 

Wave 

9775 

0.08 

— 

4 

1733-1750 

5470 

KH 

25  525 

0.16 

370  up 

5 

1751-1759 

5150 

KH 

11  275 

0.25 

850-1500 

6 

1801-1816 

4840 

KH 

22  525 

0.25 

2800-3500 

7 

1818-1825 

4620 

KH 

11  525 

0.42 

500-1000 

8 

1827-1833 

4240 

Turb 

9025 

0.39 

— 

9 

1840-1845 

3930 

Turb 

6775 

0.38 

380-550 

10 

1846-1854 

3620 

Turb 

11  275 

0.60 

490  up 

11 

1858-1903 

3300 

Turb 

7025 

0.69 

950  up 

12 

1905-1909 

3000 

Turb 

6775 

— 

860  up 

13 

1914-1918 

2670 

Turb 

5525 

— 

— 

14 

1920-1924 

2390 

Turb 

7025 

— 

— 

15 

1929-1932 

2080 

Turb 

5525 

— 

150-400 

16 

1933-1936 

1790 

Turb 

5275 

— 

— 

15  km  (i.e.,  in- valley  legs)  to  80  km  (above-mountaintop 
legs).  The  average  true  airspeed  of  UWKA  is  approxi¬ 
mately  100-110  m  s_1 for  the  long  legs  above  the  moun- 
taintop  and  approximately  80-90  m  s-1  for  the  short 
in-valley  legs.  The  data  frequency  is  25  Hz.  In  addition  to 
the  UWKA  data,  radiosondes  launched  from  the  National 
Center  for  Atmospheric  Research  (NCAR)  Integrated 
Sounding  System  (ISS)  located  at  the  Independence 
Airport  are  analyzed  as  well. 

In  this  study,  we  focus  on  the  14  March  2006  morning 
flight  executed  during  IOP  4  (referred  to  as  IOP  4a 
hereafter).  UWKA  took  off  at  approximately  1600  UTC 
[i.e.,  0800  Pacific  Standard  Time  (PST)]  and  returned  to 
the  Bishop  airport  at  approximately  2000  UTC  (i.e., 
1200  PST).  The  UWKA  sampled  from  the  middle  tro¬ 
posphere  down  to  approximately  600  m  above  the  valley 
floor  over  a  2.5-h  period  and  the  vertical  spacing  be¬ 
tween  two  adjacent  legs  was  approximately  300  m  with 
the  exception  of  the  top  two  legs.  In  addition,  we  have 
analyzed  the  morning  flight  on  2  March  during  IOP  1 
(hereafter  referred  to  as  IOP  la)  and  the  results  are 
compared  with  IOP  4a  in  section  6.  The  UWKA  in¬ 
strumental  error  ranges  are  0.5  K  for  potential  temper¬ 
ature,  1  m  s-1  for  horizontal  velocity,  0.15  m  s-1  for 
vertical  velocity,  and  0.5  hPa  for  pressure. 

b.  Mean  flow  conditions 

Results  from  the  Coupled  Ocean-Atmospheric  Me- 
soscale  Prediction  System  (CO AMPS)1  simulations  with 


1  CO  AMPS  is  a  registered  trademark  of  the  Naval  Research 
Laboratory. 


Fig.  2.  Profiles  of  (left)-(right)  U,  V,  6 ,  N2,  and  /  are  shown.  The 
dashed  profiles  are  derived  from  the  1800  UTC  14  Mar  2006 
sounding  launched  from  the  MAPR  site  in  Owens  Valley.  The 
crosses  represent  data  points  from  the  means  of  each  UWKA  leg, 
and  the  solid  curves  correspond  to  sixth-order  regressions  of  U,  V, 
and  6  using  the  UWKA  means. 

a  27-km  grid  increment  (not  shown)  indicate  that  the 
large-scale  flow  conditions  for  IOPs  la  and  4a  are  both 
characterized  by  a  pressure  trough  positioned  off  the 
Alaska  Coast,  which  directs  the  southwesterly  flow  to¬ 
ward  the  Sierra  Nevada  and  Owens  Valley.  During  the 
observational  periods,  the  prevailing  winds  at  the  moun- 
taintop  level  are  southwesterly  and  the  UWKA  flight 
tracks  are  oriented  approximately  along  the  mountaintop 
level  wind  direction.  In  the  remainder  of  the  paper,  U  and 
V  denote  the  along-track  (also  referred  to  as  longitudinal 
or  cross-valley)  and  the  normal  (also  referred  to  as  along- 
valley)  wind  components,  respectively. 

Figure  2  shows  profiles  derived  from  the  1800  UTC 
14  March  sounding  launched  from  the  NCAR  ISS  near 
Independence;  the  corresponding  variables  averaged  from 
each  leg  of  the  UWKA  IOP  4a  flight  are  included  for 
comparison.  The  radiosonde  measurements  have  been 
averaged  over  a  vertical  distance  of  300  m,  which  is 
comparable  to  the  average  vertical  distance  between  two 
adjacent  UWKA  legs.  The  wind  and  potential  tempera¬ 
ture  profiles  from  the  two  different  datasets  agree  with 
each  other  reasonably  well.  The  discrepancies  could  be 
due  to  the  horizontal  variation  of  the  measured  fields 
and  the  time  lag  between  the  radiosonde  and  the  UWKA 
flight.  For  example,  the  potential  temperature  measured 
by  the  UWKA  at  the  1790-m  level  near  1930  UTC  is 
approximately  8  K  warmer  than  that  derived  from  the 
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1800  UTC  (i.e.,  launched  near  1700  UTC)  radiosonde, 
probably  because  of  the  warming  of  the  low-level  valley 
air  associated  with  the  surface  heating  between  1700  and 
1930  UTC.  The  cross- valley  wind  profile  is  characterized 
by  weak  winds  (<10  m  s-1)  below  3.5  km  MSL,  strong 
winds  (—40  m  s-1)  above  5.5  km,  and  a  shear  layer  lo¬ 
cated  between  3.5  and  5.5  km.  The  along-valley  wind 
component  is  southerly  in  the  valley  and  is  much  weaker 
than  the  along-track  component  above  the  mountaintop 
level.  A  stable  layer  exists  between  4  and  5.5  km,  co¬ 
incident  in  altitude  with  the  shear  layer  above  the  moun¬ 
taintop.  The  air  in  the  valley  is  much  less  stable. 

Also  shown  in  Fig.  2  are  the  buoyancy  frequency 
squared  N2(z)  defined  by  N 2  =  ( g/6)(d0/dz ),  and  the  Scorer 
parameter  l(z )  defined  by  l2  =  N2/U 2  —  (1  IU)(d2Uldz2), 
where  g  is  the  gravitational  acceleration  and  0  is  the  po¬ 
tential  temperature.  It  is  noteworthy  that  N 2  computed 
from  UWKA  leg-average  potential  temperature  is  neg¬ 
ative  between  2.5  and  3  km  MSL,  indicative  of  the  pres¬ 
ence  of  an  elevated  mixed  layer.  A  similar  elevated  mixed 
layer  was  observed  during  the  Sierra  Rotor  Project  (Jiang 
and  Doyle  2008)  associated  with  a  shallow  layer  of  west¬ 
erly  flow  across  the  valley  floor.  The  N 2  profiles  derived 
from  both  the  radiosonde  measurements  and  UWKA  leg 
means  exhibit  a  maximum  between  4.5  and  5.5  km.  It  is 
also  interesting  to  note  that  the  value  of  N 2  derived  from 
the  radiosonde  is  larger  than  3  X  10-4  s-1  between  4.75 
and  5.25  km,  coincident  with  a  vertical  wind  shear  max¬ 
imum  of  approximately  0.015  s-1.  Associated  with  the 
stability  variation  and  the  increase  in  U ,  the  Scorer  pa¬ 
rameter  decreases  gradually  with  altitude  in  the  4-7-km 
layer,  which  may  serve  as  a  wave  duct  and  promote  the 
formation  of  trapped  waves  in  the  lee  of  topography 
(Scorer  1949). 

3.  Waves,  KH  instability,  and  turbulence 

The  25-Hz  vertical  velocities  observed  along  verti¬ 
cally  stacked  legs  as  a  function  of  distance  from  a  ref¬ 
erence  point  over  the  Sierra  Ridge  (i.e.,  the  western  end 
of  the  longest  leg)  are  shown  in  Fig.  3.  Within  the  vertical 
cross  section,  UWKA  sampled  small-scale  turbulence, 
kilometer-scale  eddies,  and  mountain  waves  with  wave¬ 
lengths  on  the  scale  of  20-30  km.  In  this  study,  we  use 
the  terms  “waves,”  “large  eddies,”  and  “turbulence”  to 
loosely  refer  to  perturbations  with  scales  greater  than 
6  km,  between  500  m  and  6  km,  and  smaller  than  500  m, 
respectively.  The  choice  of  the  eddy  size  range  is  largely 
based  on  the  spectral  and  wavelet  analysis  in  the  fol¬ 
lowing  sections.  For  the  convenience  of  discussion,  we 
separate  the  16  legs  into  three  groups.  The  top  three  legs 
above  5.5  km,  characterized  by  relatively  smooth  waves, 
are  referred  to  as  wave  legs.  The  flight  segments  below 
the  mountaintop  are  predominately  turbulence  and  are 


Fig.  3.  The  25-Hz  vertical  velocity  is  plotted  as  a  function  of  the 
horizontal  distance  from  a  reference  point  over  the  Sierra  ridge 
(the  west  end  point  of  the  top  leg)  for  legs  (top)-(bottom)  1-10  and 
16,  including  topography.  The  corresponding  vertical  displacement 
and  filtered  vertical  velocity  curves  are  shown  as  dashed  and  gray 
curves,  respectively.  Note  the  variation  in  ranges  of  w  and  17. 


therefore  referred  to  as  turbulence  legs.  The  four  legs 
between  4.5  and  5.5  km  MSL  are  characterized  by  coex¬ 
istence  of  smooth  waves,  large  eddies,  and  high-frequency 
turbulence.  We  refer  to  these  legs  as  KH  legs  based  on 
the  discussion  presented  later  in  the  paper. 

a.  Trapped  waves 

Also  included  in  Fig.  3  are  the  vertical  velocity  smoothed 
with  a  250-point  (i.e.,  —1  km)  running  mean  filter  for  each 
leg,  and  the  vertical  displacement  for  the  top  10  legs.  The 
vertical  displacement  r](x)  is  derived  using  the  steady- 
state  linear  relation  (Smith  et  al.  2008) 

r](x)  =  f  (w/U)dx,  (1) 

Jo 

where  w  is  the  vertical  velocity,  U  is  based  on  a  five-point 
average,  and  r](x)  is  computed  only  for  legs  with  the 
minimum  U  greater  than  3  m  s-1.  Note  that  the  curves 
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for  the  three  wave  legs  above  5.5  km  shown  in  the  top 
panel  agree  with  each  other  well. 

Regardless  of  the  time  lag  between  these  legs,  verti¬ 
cally  coherent  wave  patterns  are  evident  for  the  wave, 
KH,  and  top  turbulence  legs  near  the  mountaintop  level, 
indicating  that  the  waves  are  relatively  steady  over  the 
approximately  2.5-h  observational  period  and  penetrate 
into  the  valley.  At  least  two  wave  crests  are  evident  from 
the  wave  and  KH  legs  with  a  wavelength  of  approxi¬ 
mately  25  km  (Fig.  3).  The  first  wave  crest  is  located  over 
the  lee  slope  of  the  Sierra  ridge  and  the  second,  with 
a  much  reduced  amplitude,  is  located  over  the  eastern 
portion  of  the  valley.  The  observed  wave  crests  and  troughs 
are  vertically  aligned  (Fig.  3),  suggesting  that  these  waves 
are  probably  trapped  waves.  As  demonstrated  in  Smith 
et  al.  (2008),  a  useful  parameter  for  diagnosing  trapped 
waves  from  aircraft  data  is  the  wave  energy  equipartition 
ratio JEQR);  EQR  =  PE/KE,  where  PE  =  -(pg/2T)oy 
and  KE  =  (p/2)(w'2  +  v'2  +  w'2)  are  the  transect-averaged 
potential  and  kinetic  energy,  p  and  T  are  the  air  density 
and  temperature,  respectively,  the  overbars  denote  tran¬ 
sect  means,  and  the  primed  variables  represent  pertur¬ 
bations  relative  to  their  leg  means.  For  a  pure  propagating 
wave,  the  wave  EQR  should  be  unity  and  any  departure 
from  unity  indicates  the  presence  of  both  up-  and  down¬ 
going  waves  (Smith  et  al.  2008).  The  EQR  values  esti¬ 
mated  from  the  top  three  legs  are  0.067,  0.11,  and  0.48, 
respectively,  further  substantiating  that  the  observed 
waves  are  indeed  trapped.  The  decay  of  trapped  waves 
downstream  could  be  caused  by  boundary  layer  ab¬ 
sorption  (Jiang  et  al.  2006)  or  wave  energy  leaking  aloft. 

b.  KH  instability 

According  to  linear  theory,  KH  instability  may  occur  in 
stably  stratified  shear  flow  when  the  gradient  Richardson 
is  less  than  Va  (Miles  1961;  Howard  1961).  The  develop¬ 
ment  of  KH  instability  usually  leads  to  the  formation  of 
a  series  of  parallel  roll  vortices  (often  referred  to  as  KH 
waves,  KH  billows,  or  KH  eddies)  oriented  along  the  wind 
direction.  KH  instability  has  been  the  subject  of  numerous 
theoretical,  experimental  (e.g.,  Thorpe  1968),  and  nu¬ 
merical  (e.g.,  Klaassen  and  Peltier  1985;  Fritts  et  al.  1996) 
studies.  KH  eddies  can  be  made  visible  by  clouds  (i.e., 
billow  clouds)  or  can  be  observed  by  Doppler  radars  in 
the  optically  clear  atmosphere.  However,  direct  obser¬ 
vations  of  KH  eddies  by  research  aircraft  are  still  rare. 

Legs  4-7  are  located  between  4.5  and  5.5  km,  a  layer 
characterized  by  strong  vertical  wind  shear  (i.e.,  mean 
shear  —0.013  s-1)  and  relatively  large  static  stability 
(N2  >  2  X  10-4  s-2;  Fig.  2).  We  estimate  the  gradient 
Richardson  number  Ri,  defined  as  Ri  =  N2(UZ  +  Lz)_1, 
using  sixth-order  polynomial  regressions  of  the  transect 
means,  (i.e.,  U ,  V,  and  6)  to  obtain  the  vertical  derivatives 


Fig.  4.  Profiles  of  (left)-(right)  gradient  and  flux  Richardson 
numbers,  leg-averaged  vertical  velocity  (m  s_1),  sensible  heat  flux 
(W  m-2),  vertical  flux  of  the  cross- valley  momentum  component 
(N  m~2),  and  TKE  (m2  s-2)  for  each  leg. 

of  those  variables  (Fig.  4).  In  addition,  the  flux  Richardson 
number  Ri f  =  bis ,  where  b  =  (g/6)0"w"  is  the  buoyancy 
production  rate  of  turbulent  kinetic  energy  (TKE)  and 
s  =  u"w"(dU/dz )  +  v"w"(dV/dz)  is  the  shear  production 
rate,  estimated  from  all  flight  legs  is  included  in  Fig.  4  as 
well.  The  double-primed  variables  denote  perturbations 
associated  with  large  eddies  and  turbulence,  obtained 
by  subtracting  the  6-km  moving  averages  to  remove 
mountain  wave  signals.  Both  the  gradient  and  the  flux 
Richardson  numbers  are  less  than  unity  below  5  km  MSL 
and  increase  to  2  near  5.5  km  MSL.  At  the  levels  of  the 
four  KH  legs,  the  mean  Richardson  numbers  are  rela¬ 
tively  small  but  still  larger  than  the  critical  Richardson 
number  (i.e.,  1/a).  This  is  consistent  with  the  nearly  lami¬ 
nar  waves  observed  outside  the  first  crest  along  the  KH 
legs.  The  occurrence  of  KH  eddies  in  the  first  wave  crest 
likely  results  from  the  interaction  between  trapped  waves 
and  the  stable  shear  layer  between  4.5  and  5.5  km.  As 
demonstrated  in  the  appendix,  a  trapped  wave  can  modify 
the  local  Richardson  number  by  changing  the  vertical 
wind  shear,  static  stability  and  vertical  acceleration.  In 
general,  the  Richardson  number  tends  to  be  reduced  in 
wave  crests  and  the  reduction  is  proportional  to  the  wave 
amplitude.  Although  two  wave  crests  are  observed  along 
these  KH  legs,  the  amplitude  of  the  first  wave  is  sub¬ 
stantially  larger  than  the  second.  Consequently,  the  local 
Richardson  number  should  be  smaller  in  the  first  wave 
crest  where  KH  eddies  and  turbulence  are  observed. 

c.  Turbulence 

Although  the  in-valley  flight  legs  are  predominantly 
turbulent  in  character,  the  filtered  vertical  velocity  still 
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reveals  the  existence  of  wave  signatures  (>6  km),  es¬ 
pecially  along  the  upper  in-valley  legs  (i.e.,  above  3  km 
MSL).  Large  eddies  are  also  apparent  at  all  levels,  in¬ 
dicative  of  complex  multiscale  interaction  between  tur¬ 
bulence  and  waves. 

The  mean  vertical  velocity  w,  sensible  heat  flux 
Sflx  (: =pCp6W '),  momentum  Uflx  ( =pu"w ")  flux,  and  tur- 
bulence  kinetic  energy  per  unit  mass  e  [=(u"2  +  v"2  +  w"2)l 
2;  hereafter  referred  to  as  TKE  for  simplicity]  are  also 
shown  in  Fig.  4.  Above  the  mountaintop  level,  the  leg-mean 
vertical  velocity  is  negative  and  on  the  order  of  0.5  m  s-1, 
associated  with  mountain  wave-induced  descent  in  the 
lee  of  the  Sierra  ridge.  In  Owens  Valley,  the  wave-induced 
mean  vertical  velocity  is  positive,  probably  associated  with 
the  ascending  branch  of  the  mountain  waves. 

Above  5.5  km  MSL,  the  sensible  heat  and  cross-valley 
momentum  fluxes  (wave  contribution  has  been  re¬ 
moved),  and  TKE  are  nearly  zero,  which  is  consistent 
with  the  laminar  waves  observed  along  the  top  three  legs. 
For  the  KH  and  in- valley  turbulence  flight  legs,  the  av¬ 
erage  sensible  heat  fluxes  are  negative,  indicative  of 
downward  heat  transfer  by  turbulence.  The  largest 
downward  heat  flux  is  observed  along  the  lowest  KH  leg 
(i.e.,  —4.6  km  MSL),  where  the  air  is  stable  and  turbu¬ 
lence  is  strong.  The  cross- valley  momentum  flux  is  nega¬ 
tive  in  the  valley,  implying  downward  momentum  mixing. 
As  expected,  the  momentum  flux  amplitude  also  reaches 
a  maximum  approximately  at  4.6  km  MSL,  where  strong 
vertical  wind  shear  and  turbulence  are  present.  A  sec¬ 
ondary  maximum  is  located  below  the  ridge  crest,  prob¬ 
ably  associated  with  wave-induced  vertical  wind  shear. 
The  TKE  profile  shows  a  maximum  (—39  m2  s-2)  near 
4.6  km  MSL,  apparently  associated  with  KH  eddies.  In 
the  valley,  TKE  shows  little  variation  with  height,  which 
reflects  the  presence  of  a  deep  turbulent  layer  in  Owens 
Valley. 

4.  Breaking  KH  billows 

The  high-frequency  variations  superposed  on  the 
kilometer-scale  eddies  observed  along  legs  5-7  suggest 
that  these  eddies  are  overturning  or  breaking.  In  this 
section,  we  will  further  examine  characteristics  of  these 
large  eddies. 

Within  the  first  wave  crest,  kilometer-scale  eddies 
oriented  along  the  prevailing  wind  direction  are  evident 
in  all  four  KH  leg  measurements  (Fig.  3).  However,  the 
characteristics  of  small-scale  turbulence  observed  along 
the  four  legs  are  substantially  different.  The  top  KH 
leg  (i.e.,  leg  4)  is  characterized  by  four  vertical  velocity 
maxima  within  the  first  wave  crest  with  a  horizontal 
wavelength  of  approximately  2-3  km.  Perturbations  with 
wavelengths  less  than  500  m  are  negligible.  Along  leg  5, 


the  kilometer-scale  eddies  are  relatively  smooth  as  well 
except  for  a  few  narrow  turbulent  zones  (Fig.  3).  The 
lower  two  KH  legs  (i.e.,  6  and  7)  are  significantly  more 
turbulent  throughout  the  first  wave  crest.  The  coexistence 
of  kilometer-scale  eddies  and  small-scale  turbulence 
observed  along  legs  5-7  is  consistent  with  simulated 
breaking  KH  eddies  in  previous  studies  (e.g.,  Klaassen 
and  Peltier  1985;  Fritts  et  al.  1996).  The  difference  in  the 
characteristics  of  small-scale  turbulence  along  legs  4-7 
could  be  due  to  the  intersection  of  the  flight  segments 
with  the  different  portions  of  these  breaking  KH  eddies. 
According  to  previous  numerical  studies  (e.g.,  Klaassen 
and  Peltier  1985;  Fritts  et  al.  1996),  breaking  KH  billows 
(eddies)  are  characterized  by  turbulent  billow  cores  con¬ 
nected  by  less  turbulent  braids  where  secondary  instability 
sometimes  occurs.  If  the  aircraft  flight  level  is  above  the 
turbulent  billow  cores  (e.g.,  leg  4),  only  large  eddies  are 
observed.  For  a  flight  leg  that  cuts  through  the  top  of  the 
billow  cores  (e.g.,  leg  5),  the  aircraft  samples  large 
eddies  separated  by  narrow  turbulent  zones  where  the 
flight  leg  intersects  with  the  billow  cores.  The  lower  two 
flight  legs  (e.g.,  legs  6  and  7)  observe  KH  eddies  super¬ 
posed  with  strong  turbulence  as  they  pass  through  tur¬ 
bulent  billow  cores. 

To  examine  the  structure  and  phase  relations  of  these 
multiscale  perturbations,  the  vertical  velocity,  cross¬ 
valley  horizontal  velocity,  and  potential  temperature 
along  legs  5-7  over  a  15-km  segment  within  the  first 
wave  crest  are  shown  in  Fig.  5.  It  is  noteworthy  that 
variations  of  the  along- valley  velocity  component  (i.e., 
normal  to  the  flight-track  component;  not  shown)  asso¬ 
ciated  with  kilometer-scale  eddies  are  much  smaller 
than  the  cross-valley  component.  For  leg  5,  large  eddies 
with  a  horizontal  scale  of  approximately  2  km  are  evi¬ 
dent  in  both  the  vertical  velocity  and  the  cross-valley 
wind  component  measurements  (Fig.  5a).  The  corre¬ 
sponding  vertical  velocity  range  (i.e.,  estimated  from  the 
filtered  curve)  is  approximately  6  ms”1  and  the  corre¬ 
sponding  U  variations  are  much  smaller  (~3  m  s"1). 
There  seems  to  be  an  approximate  90°  phase  lag  be¬ 
tween  w  and  U  variations  for  these  large  eddies.  The 
potential  temperature  variation  is  nearly  180°  out  of 
phase  with  the  U  variation  and  90°  out  of  phase  with  the 
vertical  velocity.  Although  these  KH  eddies  are  highly 
nonlinear,  linear  theory  may  provide  some  useful  quali¬ 
tative  insights  into  these  phase  relations.  For  a  given 
horizontal  wavenumber  k  the  linear  polarization  re¬ 
lations  between  w,  U ,  and  0  perturbations  can  be  written 
as  (e.g.,  Nielsen  1992) 

w{k)  =  -3^[c.  +  i(U  -  cr)]0(k)  and  (2) 


July  2010 


JIANG  ET  AL. 


2361 


10  15  20  25 


Distance  (km) 

Fig.  5.  Plots  of  w,  U,  and  6  as  a  function  of  distance  for  the 
segment  between  10  and  25  km  for  KH  legs  (a)  5,  (b)  6,  and  (c)  7. 
The  corresponding  filtered  curves  are  also  included. 


u(k)  =  -jw(k),  (3) 

respectively,  where  c  =  cr  +  ict  is  the  complex  wave 
phase  speed  and  (w,  u,  6)  denote  the  corresponding 
Fourier  components  of  the  (w,  U ,  0)  variations.  The 
vertical  wavelength  m  is  given  by  m2  =  [N2/(U  —  c )2]— 
[(d2U/dz2)/(U-c)]-k2. 


Based  on  (2)  and  (3),  the  vertical  velocity  and  po¬ 
tential  temperature  variations  should  be  90°  out  of  phase 
for  a  neutral  wave  mode  (i.e.,  ct  =  0),  and  for  a  neutral 
and  propagating  (i.e.,  q  =  0  and  m 2  >  0)  wave,  the  phase 
lag  between  the  vertical  and  horizontal  velocity  pertur¬ 
bations  should  be  180°  (or  -180°).  The  phase  relations 
for  kilometer-scale  eddies  observed  along  leg  5  are 
consistent  with  a  propagating  neutral  plane  wave  with 
the  phase  line  tilted  westward  with  height. 

Along  leg  5,  sharp  perturbations  appear  in  three 
narrow  zones  (-0.2-0.5  km  wide)  located  near  17,  19, 
and  21.5  km,  respectively,  and  the  maximum  w  range  of 
the  high-frequency  perturbations  reaches  approximately 
16  m  s”1  (Fig.  5a).  In  these  turbulence  zones,  the  verti¬ 
cal  velocity  perturbations  associated  with  large  eddies  are 
nearly  zero  and  the  corresponding  6  perturbations  reach 
local  minima,  which  is  also  consistent  with  the  findings  of 
Nielsen  (1992)  from  aircraft  observations  of  KH  waves 
associated  with  a  coastal  frontal  inversion.  He  noticed 
that  velocities  are  “turbulent  during  the  cold  half  of  the 
wave  and  smoother  during  the  warm  half  of  the  wave.” 
For  leg  6,  several  large  eddies  with  two  predominant 
horizontal  scales  of  approximately  1  and  2  km  are  em¬ 
bedded  in  strong  turbulence  signals.  Similarly,  leg  7 
is  characterized  by  both  strong  turbulence  and  large 
eddies  with  a  predominant  horizontal  scale  of  approxi¬ 
mately  3.5  km.  In  general,  for  these  kilometer-scale  eddies, 
the  U  and  6  variations  appear  to  be  approximately  in  phase 
with  each  other  and  approximately  180°  out  of  phase  with 
the  w  variations.  These  phase  relations  imply  negative 
momentum  and  sensible  heat  fluxes  associated  with  these 
eddies,  the  impact  of  which  on  large  eddies  can  be  seen 
from  the  TKE  equation  for  a  two-dimensional  com¬ 
pressible  flow: 


De 

~Dt 


-  -u'w' —  +  -w'0'  - 

dz  e 


dw'e' 

dz 


1  dw'p’ 
p  dz 


(4) 


where  the  right-hand  terms  correspond  to  shear  pro¬ 
duction,  buoyancy  production,  vertical  turbulent  trans¬ 
port,  pressure  correlation,  and  dissipation  rate  of  TKE, 
respectively.  According  to  (4),  given  a  positive  shear  in 
the  mean  winds,  a  negative  momentum  flux  indicates 
a  positive  TKE  production  rate  and  a  negative  heat  flux 
tends  to  inhibit  turbulence  development.  The  turbulent 
transport  and  pressure  correlation  terms  are  typically 
small  and  difficult  to  estimate.  To  evaluate  the  relative 
importance  of  the  first  two  terms  for  large  eddies,  we 
define  R(k )  as  the  ratio  of  the  first  two  terms  in  their 
spectral  forms;  that  is, 


m= 


gCJV(d_uXl 

eCvu&Ate) 


(5) 
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Table  2.  Coh^,  Coh^,  </>wu,  4>we,  Cwu,  and  Cw6  for  A  derived  from  flight  legs  5,  6,  and  7.  The  corresponding  R(k  =  2ttIX)  values  are 

included  as  well. 


Leg 

U  (m  s-1) 

A  (km) 

Ku  O 

CohWM 

Cwu  (m2  s  x) 

4>wd  (°) 

Coh^ 

CwG  (K  m) 

m 

5 

28 

2.0 

-105 

0.98 

-17.8 

110 

0.80 

-8.1 

1.5 

6 

27 

1.3 

-170 

0.92 

-170 

-155 

0.73 

-12 

0.24 

7 

24 

3.5 

165 

0.80 

-250 

160 

0.83 

-50 

0.66 

where  Cw0  and  Cwu  denote  the  cospectra  of  vertical 
velocity,  potential  temperature,  and  cross-valley  veloc¬ 
ity  component  variations.  Note  that  R(k)  resembles  a 
Richardson  number  and  measures  the  relative  ampli¬ 
tudes  of  the  shear  and  buoyancy  production  of  pertur¬ 
bation  energy  for  a  given  k.  When  R(k )  <  1,  the  shear 
production  term  dominates  and  eddies  with  k  could  grow 
by  extracting  mechanical  energy  from  the  mean  flow  to 
overcome  the  stratification  effect.  For  legs  5-7,  the 
dominant  wavelengths  A  =  2i rlk  and  the  corresponding 
R(k)  values  are  calculated  using  Fourier  transform  and 
listed  in  Table  2.  Taylor’s  hypothesis  of  “frozen”  tur¬ 
bulence  is  used  to  convert  the  UWKA  measurements 
from  the  frequency  space  into  the  wavenumber  space.  It 
is  noteworthy  that  the  UWKA  tends  to  maintain  a  con¬ 
stant  true  airspeed  along  each  flight  leg;  as  a  result,  the 
ratio  between  the  longitudinal  velocity  variations  and 
the  UWKA  true  airspeed  is  on  the  order  of  10%  or  less. 
Therefore,  although  the  large  velocity  variations  associ¬ 
ated  with  large -amplitude  waves  and  turbulence  observed 
during  IOP  4a  may  have  some  effect  on  the  accuracy  of 
the  estimated  dominant  eddy  sizes,  the  Taylor  hypothesis 
is  still  valid.  The  dominant  wavenumbers  are  selected  by 
identifying  energy-containing  peaks  in  velocity  and  po¬ 
tential  temperature  power  spectra  using  the  25-Hz  air¬ 
craft  data  over  the  10-25 -km  segment  in  the  first  wave 
crest  after  removing  the  linear  trend.  It  is  interesting  that 
the  U,  w,  and  6  variations  are  highly  coherent  for  these 
wavenumbers.  The  corresponding  V  (i.e.,  along-valley 
component)  variations  are  much  smaller  and  the  co¬ 
herence  (Coh)  between  V  and  w  variations  is  generally 
below  0.4.  This  is  consistent  with  two-dimensional  KH 
billows,  the  axis  of  which  is  perpendicular  to  the  wind 
shear. 

For  all  three  legs,  the  buoyancy  tends  to  inhibit  the 
growth  of  the  identified  eddies,  as  expected  in  a  stably 
stratified  flow.  The  phase  relations  between  U  and  w  al¬ 
low  eddies  to  extract  energy  from  the  positively  sheared 
flow.  For  leg  5,  the  buoyancy  effect  dominates  and  the 
identified  eddies  are  unable  to  grow  even  if  the  dissipa¬ 
tion  is  negligible.  Again,  this  is  consistent  with  the  hy¬ 
pothesis  that  leg  5  passes  through  the  top  of  breaking  KH 
eddies.  Along  legs  6  and  7  the  R(k)  values  are  less  than 
unity,  indicating  that  the  shear  production  of  energy  is 
large  enough  to  overcome  the  buoyancy  effect  and  the 


eddies  are  probably  growing.  It  is  also  interesting  that, 
using  linear  relations  (2)  and  (3)  and  the  phase  angles 
between  w  and  6  in  Table  2,  we  obtain  that  ct  >  0  for 
both  legs  6  and  7,  that  U  -  cr>  0  (i.e.,  disturbance  drifts 
with  the  local  flow)  for  leg  6,  and  that  U  —  cr<  0  for  leg 
7,  implying  that  the  steering  level  is  probably  located 
between  legs  6  and  7.  However,  we  should  be  cautious 
about  applying  linear  theory  to  the  observed  KH  eddies, 
which  are  highly  nonlinear,  especially  for  legs  6  and  7. 
Finally,  we  note  that  the  observed  largest  eddy  size  is 
3.5  km  and  the  depth  of  the  stable  layer  is  approximately 
0.5  km,  and  the  ratio  of  the  two  is  comparable  to  that  of 
KH  waves  from  previous  studies  (i.e.,  —7.5;  Miles  and 
Howard  1964). 

In  summary,  the  characteristics  of  the  kilometer-scale 
eddies  observed  along  legs  4-7  are  consistent  with 
breaking  KH  billows  described  in  previous  studies  (e.g., 
Klaassen  and  Peltier  1985;  Nielsen  1992;  Fritts  et  al.  1996). 
These  eddies  are  not  likely  generated  by  breaking  trapped 
waves.  The  wave  nonlinearity  index  for  trapped  waves — 
a  =  2V| T7|max/T7,  where  |i7|max  is  the  maximum  vertical 
displacement  for  each  leg — is  relatively  small  for  the  wave 
and  KH  legs  (Table  1),  implying  that  nonlinear  processes 
such  as  wave  steepening  and  breaking  contribute  little  to 
the  eddy  and  turbulence  generation. 

5.  Characteristics  of  large  eddies  and  turbulence 

The  preliminary  analysis  of  the  KH  and  turbulence 
legs  in  the  previous  section  reveals  the  coexistence  of 
waves,  large  eddies,  and  small-scale  turbulence  from  the 
lower  to  midtroposphere  over  Owens  Valley.  In  this  sec¬ 
tion,  we  examine  further  the  characteristics  of  the  observed 
eddies  and  turbulence  using  three  different  methods, 
namely,  Fourier  transform-based  spectral  analysis,  wave¬ 
let  analysis,  and  structure  functions. 

a.  Spectral  analysis 

Previous  studies  have  demonstrated  that  Fourier 
spectral  analysis  is  a  useful  tool  for  scale  separation  and 
for  identifying  turbulence  inertial  subranges  (e.g.,  Lilly 
and  Lester  1974;  Nastrom  and  Gage  1985,  among  others). 
As  an  example,  the  power  spectra  of  the  vertical  and 
longitudinal  velocity  components  and  the  potential  tem¬ 
perature  for  leg  5  are  shown  in  Fig.  6a.  The  25-Hz  data 


July  2010 


JIANG  ET  AL. 


2363 


Frequency  (Hz) 


Fig.  6.  (a)  The  power  spectra  density  (PSD)  of  w,  U'  (m2  s_1), 
and  6  (K2  s)  as  a  function  of  frequency  in  log-coordinates  for  KH 
leg  5.  The  periods  (s)  and  horizontal  scales  (km)  corresponding  to 
three  energy-containing  peaks  are  labeled,  (b)  The  PSD  of  w  for 
the  four  KH  legs.  The  bold  lines  with  slopes  of  — 5/3  are  included  for 
comparison. 

have  been  linearly  detrended  before  applying  the  Fourier 
transform,  and  the  obtained  spectra  have  been  averaged 
over  each  equal -log  interval  [i.e.,  Alog(/)  =  0.1]  of  fre¬ 
quency.  Only  the  portion  between  0.01  and  10  Hz  is 
shown.  A  number  of  maxima  are  evident  in  the  power 
spectra  of  the  u ,  w,  and  6  variations,  implying  the  pres¬ 
ence  of  multiscale  waves  and  large  eddies  along  the  KH 
legs.  For  example,  in  the  w  spectrum  of  leg  5  there  are 


Fig.  7.  As  in  Fig.  6a,  but  for  flight  leg  11. 

three  peaks,  centered  approximately  at  8.2,  2.0,  and 
0.92  km,  respectively,  that  satisfy  the  95%  confidence  test 
(Bendat  and  Piersol  1986).  Particularly,  the  2.0-km  maxi¬ 
mum  appears  in  all  three  spectra  and  the  scale  is  compa¬ 
rable  to  the  visual  estimation  of  the  KH  eddy  wavelength 
within  the  first  wave  crest  (Fig.  3).  The  mean  spectral 
slopes  are  approximately  -5/3  over  the  high-frequency 
range,  from  10  Hz  up  to  approximately  0.05  Hz,  which 
corresponds  to  approximately  2  km  in  length  scale. 

The  w  spectra  of  the  four  KH  legs  share  some  common 
characteristics  such  as  multiple  kilometer-scale  energy- 
containing  maxima  (Fig.  6b).  For  the  three  lower  KH 
legs,  the  mean  power  spectral  slopes  are  approximately 
-5/3  for  scales  less  than  2  km,  suggestive  of  an  inertial 
subrange.  The  spectrum  of  the  top  KH  leg  (i.e.,  leg  4) 
exhibits  a  mean  slope  gentler  than  -5/3  for  higher  fre¬ 
quencies.  The  lack  of  a  well-defined  inertial  subrange  in 
the  spectra  of  the  top  KH  leg  is  consistent  with  the  weak 
high-frequency  perturbations  along  this  leg  (Fig.  3).  For 
turbulent  legs  below  the  mountaintop  level,  the  power 
spectra  follow  the  -5/3  power  law  from  10  Hz  up  to  ap¬ 
proximately  0.1  Hz  (~1  km;  Fig.  7),  and  multiple  energy- 
containing  maxima  are  evident  over  lower  frequencies, 
corresponding  to  kilometer-scale  eddies.  The  u  and  w 
spectra  follow  each  other  more  closely,  suggesting  that 
the  turbulence  is  more  isotropic  below  the  mountaintop 
level  than  at  the  KH  leg  levels. 

b.  Wavelet  analysis 

Wavelet  analysis  has  been  widely  used  for  the  study  of 
intermittent  turbulence  and  wave-turbulence  interactions 
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(e.g.,  Argoul  et  al.  1989;  Farge  1994;  Demoz  et  al.  1998). 
Using  wavelet  analysis,  localized  turbulence  structures 
can  be  readily  revealed  in  a  time  period  [(frequency)-1] 
domain.  In  the  following,  orthonormal  wavelet  trans¬ 
formations  are  applied  to  the  UWKA  data  to  examine 
the  turbulence  energy  distribution  in  both  the  flight  time 
(space)  and  frequency  for  each  leg.  The  mother  wavelet 
we  use  here  is  the  Morlet  wavelet  (Morlet  et  al.  1982), 
one  of  the  most  commonly  applied  in  turbulence  studies. 
Shown  in  Fig.  8  are  the  wavelet  power  spectra  of  the 
normalized  vertical  velocity  perturbations  w'lcrw  for  five 
KH  legs,  where  aw  is  the  standard  deviation  of  w.  We 
only  discuss  signals  above  the  white  curves,  which  cor¬ 
respond  to  the  95  %  significance  level. 

Along  these  KH  legs,  the  perturbation  energy  is  largely 
contained  in  wave-scale  motion,  associated  with  the  trap¬ 
ped  waves  (i.e.,  warm  colors  with  time  scale  longer  than 
60  s  in  Fig.  8).  For  the  top  three  KH  legs,  small-scale 
turbulence  appears  in  the  first  wave  crest,  corresponding 
to  the  breaking  KH  eddies  in  Fig.  3.  For  example,  along 
leg  5  there  are  energy  maxima  centered  near  60  s  (~8  km) 
and  16  s  (~2  km)  in  frequency,  which  is  consistent  with 
the  power  spectrum  of  w  in  Fig.  6a.  The  2-km  maximum 
splits  into  four  vertically  oriented  conelike  patches  point¬ 
ing  toward  higher  frequencies,  suggesting  downscale  en¬ 
ergy  cascade.  For  the  lower  KH  legs,  nearly  equally  spaced 
cones  are  evident  and  the  small-scale  turbulence  becomes 
progressively  stronger.  It  is  noteworthy  that  these  conelike 
energy  maxima  are  characterized  by  distinctive  maxima 
in  frequency  rather  than  by  a  gradual  decrease  toward 
higher  frequency,  which  is  consistent  with  the  presence 
of  multiple  energy-containing  peaks  in  the  power  spec¬ 
tra  (Fig.  6b).  This  implies  that  frequency  preferences  exist 
in  the  downscale  energy  cascade  associated  with  the 
breakdown  of  KH  eddies.  Figure  8  also  indicates  that  the 
turbulence  tends  to  develop  in  the  second  wave  crest 
along  the  two  lower  KH  legs.  However,  the  turbulence  is 
much  less  pronounced  in  strength  and  less  organized 
than  in  the  first  wave  crest. 

For  turbulence  legs,  small-scale  turbulence  is  nearly 
uniformly  distributed,  implying  that  the  turbulence  is 
more  fully  developed  in  the  valley  than  at  the  KH  levels. 
The  vertically  oriented  coherent  conelike  structures 
are  not  as  evident  as  in  the  KH  legs.  However,  energy- 
containing  large  eddies  (i.e.,  scales  between  0.5-6  km) 
are  evident  along  all  the  legs.  Specifically,  some  verti¬ 
cally  coherent  large  eddies  (Fig.  8)  are  located  over  the 
eastern  side  of  the  valley,  probably  associated  with  wave 
breaking. 

c.  Structure  functions 

In  this  section,  second-  and  third-order  structure  func¬ 
tions  for  each  leg  are  evaluated  to  further  examine  the 
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Fig.  8.  The  wavelet  spectra  of  w  (normalized  by  its  variance)  for 
the  four  KH  flight  legs  and  a  turbulence  leg  are  shown.  The  time  (s) 
is  labeled  from  west  to  east  on  abscissa  and  the  period  [i.e., 
(frequency)-1]  is  shown  in  a  logarithmic  scale.  The  power 
spectra  are  color  filled  as  2n,  with  n  =  0,  1,  2,  3,  4,  5,  6,  and  7, 
corresponding  to  cold  to  warm  colors.  The  white  curves  corre¬ 
spond  to  95%  significance  level.  The  flight  leg  numbers  are  la¬ 
beled.  The  locations  of  the  first  wave  crest  are  highlighted  by 
dashed  ellipses  for  the  KH  legs. 

turbulence  characteristics.  The  one -dimensional  second- 
order  structure  function  of  a  generic  variable  <fi  is  de¬ 
fined  as 

D4Jr)  =  <(S</>')2},  (6) 

where  8<fi'  =  f>'(x  +  r)  -  4>'(x),  <fif  =  <fi  -  <fi  is  the  per¬ 
turbation  of  f>  relative  to  the  mean  value  <fi,  x  is  the 
longitudinal  axis,  r  is  the  separation  distance,  and  the 
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Fig.  9.  The  normalized  second- order  structure  functions  (a) 
Dwwr~m  and  (b)  Duur~ 2/3  for  the  four  KH  legs  and  leg  11  of  IOP  4a 
are  plotted  as  a  function  of  the  separation  distance. 


angle  brackets  represent  arithmetic  averaging.  In  the 
inertial  subrange,  the  second-order  structure  functions 
of  the  velocity  components  u ,  v ,  and  w  can  be  scaled  as 
D  ,D  ,D  ~  r 2/3 .  The  second-order  structure  func- 
tions  of  the  longitudinal  and  vertical  velocity  compo¬ 
nents  normalized  by  r273  (i.e.,  Duu  =  Duur~ 2/3  and  Dww  = 
Dr~213)  are  shown  in  Fig.  9.  It  is  evident  that  both  Dutut 
and  Duu  are  approximately  constant  from  the  smallest 
resolvable  scale  up  to  1  km  for  the  lower  three  KH  legs 
and  the  turbulence  leg,  which  is  consistent  with  the 
scaling  in  the  inertial  subrange.  For  r  >  1  km,  these 
curves  drop  off  with  further  increase  of  r.  This  is  in 
general  agreement  with  the  inertial  subrange  identified 
from  the  power  spectra  in  Figs.  6  and  7.  For  the  wave  legs 
(not  shown)  and  the  top  KH  leg,  Dww(r)r~213  decreases 
with  an  increase  in  r  for  r  <  200  m,  implying  that  the 
inertial  subrange  for  the  wave  legs  is  not  resolvable  (i.e., 
less  than  10  m),  which  is  consistent  with  the  absence  of 
the  -2/3  power  law  in  the  corresponding  power  spectrum 
(Fig.  6b).  The  kilometer  scale  energy-containing  peaks 
shown  in  the  power  spectra  are  much  less  obvious  in  the 
corresponding  Dww  and  Duu.  As  pointed  out  by  Mahrt 
and  Gamage  (1987),  an  energy  peak  in  power  spectra 
tends  to  spread  over  a  range  of  separation  distances  in 
a  structure  function. 
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Fig.  10.  (a)  VAR  defined  by  (6)  and  (b)  the  buoyancy  length 
[(7)]  for  five  IOP  4a  legs  (i.e.,  4-7  and  11)  are  plotted  as  a  function 
of  separation  distance. 


A  pair  of  useful  turbulence  indices  can  be  derived 
using  the  second-order  structure  functions,  namely  the 
velocity  aspect  ratio  VAR(r)  and  the  turbulence  buoy¬ 
ancy  length  L(r),  defined  as  (Mahrt  and  Gamage  1987) 


VAR(r)  =  2Dww(r)l[Duu(r)  +  Dvv(r)\  and  (7) 


W  =  Dww 


(8) 


respectively,  where  D w  and  Dee  are  the  second-order 
structure  functions  of  the  along-valley  velocity  compo¬ 
nent  and  potential  temperature.  For  legs  7  and  11,  VAR 
is  near  unity  for  r  <  100  m  and  decreases  to  approxi¬ 
mately  0.5  for  r  =  1  km,  implying  that  the  small-scale 
turbulence  is  nearly  isotropic,  and  large  eddies  are  more 
horizontal  because  of  the  stratification  effect  (Fig.  10). 
The  velocity  aspect  ratios  are  much  smaller  for  the  wave 
legs  (not  shown)  and  increase  with  decreasing  altitude 
for  the  four  KH  legs.  The  buoyancy  length  for  the  top 
KH  leg  is  much  smaller  than  other  legs  because  of  the 
stronger  stratification  aloft  (Fig.  2).  In  general,  the  buoy¬ 
ancy  length  increases  with  decreasing  altitude  for  legs 
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Fig.  11.  The  normalized  third-order  structure  functions  are 
plotted  vs  r  for  legs  (a)  6  and  (b)  11.  The  bold  solid  curve,  red 
circles,  and  blue  crosses  correspond  to  S 3/r,  Duu  u  ,  and  —Duu  u 
(when  Duu  u  <  0),  respectively. 


1-10,  reaches  a  maximum  at  leg  11  (—3300  m),  and  then 
decreases  further  with  decreasing  altitude. 

In  the  remainder  of  this  section,  we  evaluate  three  third- 
order  structure  functions  for  each  flight  segment,  namely, 
the  third-order  absolute  structure  function  of  the  longi¬ 
tudinal  velocity  S3,  defined  as 

s3  =  <|SM(r)|3};  (9) 

the  third-order  structure  function  Duu  u  ,  defined  by 

DUUtUt  =  ((Su)3)  +  2(Su(Sut)2),  (10) 

where  (( 8uT )2)  =  [((5u)2)  +  ((5w)2)]/2;  and  the  struc¬ 
ture  skewness,  defined  as 

S  (r)  =  D  /(. D  )3/2,  (11) 

where  Dwww  =  ((8w)3).  In  the  turbulence  inertial  sub¬ 
range,  S3  is  proportional  to  r  (Cho  et  al.  2001)  and 
Duutut  =  -4(e)r/3  (Lindborg  1996;  Antonia  et  al.  1997). 
The  sign  of  DUUtUt  indicates  the  direction  of  the  turbulence 
energy  cascade  (Frisch  1995;  Lindborg  1996);  the  energy 
cascade  is  downscale  if  DUUtUt  is  negative  (i.e.,  (e)  >  0  or 
a  positive  dissipation  rate)  and  upscale  if  it  is  positive  (i.e., 
a  negative  dissipation  rate).  The  structure  skewness  is 


r  (m) 

Fig.  12.  The  structure  skewness  is  plotted  as  a  function  of  separation 
distance  for  the  four  KH  legs  and  the  turbulence  leg  11. 


negative  if  the  shear  generation  of  turbulence  dominates 
(Mahrt  and  Gamage  1987). 

Figure  11  shows  S3/r  and  DUUtUt  =  -3DUUtUt/4v  for 
the  four  KH  and  one  turbulence  flight  segments.  Over 
turbulence  inertial  subranges,  these  functions  should  be 
constant  with  r  and  Duu  u  should  be  equal  to  the  mean 
turbulence  dissipation  rate.  For  leg  6,  the  normalized 
third-order  structure  functions  are  approximately  con¬ 
stant  for  r  between  100  and  1000  m,  an  inertial  subrange 
identified  from  spectral  and  second-order  structure 
function  analysis.  It  is  noteworthy  that  the  turbulence 
dissipation  rate  is  negative  between  2900  and  3500  m 
(Fig.  11a),  indicative  of  a  possible  upscale  energy  transfer. 
The  negative  dissipation  rates  are  found  for  some  of  the 
turbulence  legs  as  well  (Table  1).  For  example,  for  leg  11, 
the  turbulence  dissipation  rate  is  negative  for  scales  of 
1  km  or  longer  (Fig.  lib). 

The  KH  legs  are  characterized  by  large  negative 
structure  skewness  (Fig.  12),  suggesting  that  the  turbu¬ 
lence  above  the  mountaintop  is  primarily  shear  generated, 
which  is  consistent  with  the  velocity  and  Richardson 
number  profiles  (Figs.  2  and  4).  From  the  smallest  re¬ 
solvable  scale  to  1  km,  the  skewness  functions  for  the 
upper-turbulence  legs  are  negative  as  well,  with  much 
smaller  amplitude  than  the  KH  legs,  associated  with 
a  weaker  but  well-defined  forward  wind  shear.  For  some 
low-turbulence  legs  (—2500  m  or  lower;  not  shown)  the 
skewness  function  oscillates  near  zero,  indicating  that  the 
shear  effect  is  weak  and  buoyancy  generation  of  turbu¬ 
lence  may  be  important  at  these  levels. 

Finally,  it  is  worth  noting  that  turbulence  over  complex 
terrain  is  inhomogeneous  in  nature.  Consequently,  the 
average  over  each  flight  leg,  with  substantial  variability  in 
turbulence  characteristics,  may  not  necessarily  represent 
local  turbulence  statistics  in  a  quantitative  sense.  This  is 
especially  true  for  the  spectral  and  structure  function 
analysis  and  is  less  of  a  problem  for  wavelet  analysis. 
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6.  Comparisons  between  IOPs  la  and  4a 

Large  eddies  embedded  in  shear  layers  and  trapped 
waves  have  been  documented  during  several  other 
T-REX  flights  in  addition  to  IOP  4a.  Among  them,  IOP 
la  is  the  only  flight  that  documented  vertically  coherent 
KH  eddies  along  three  or  more  flight  legs  above  the 
mountaintop  level.  For  the  sake  of  generality,  in  this 
section  we  will  examine  the  IOP  la  flight  and  discuss 
similarities  and  differences  of  the  turbulence  characteris¬ 
tics  observed  during  these  two  flights. 

The  large-scale  flow  conditions  for  IOPs  la  and  4a  are 
quite  similar.  The  cross-valley  wind  profile  is  character¬ 
ized  by  weak  winds  below  3  km  MSL,  a  strong  shear  layer 
between  3  and  5  km,  and  strong  winds  aloft  (Fig.  13). 
In  the  shear  layer,  the  mean  and  maximum  vertical  wind 
shears  are  0.01  and  0.015  s-1,  respectively,  and  the  gradi¬ 
ent  Richardson  number  is  less  than  unity.  The  atmosphere 
has  greater  static  stability  below  3  km  and  becomes 
gradually  less  stable  aloft.  Consequently,  the  Scorer  pa¬ 
rameter  exhibits  a  sudden  decrease  between  3  and  4  km, 
which  favors  the  formation  of  trapped  waves. 

UWKA  observed  relatively  steady  trapped  waves 
as  indicated  by  the  vertically  coherent  wave  patterns 
(Fig.  14).  The  EQR  values  computed  from  wave  legs  1-4 
are  respectively  0.99,  0.54,  0.37,  and  0.32,  implying  that 
the  waves  are  trapped  below  approximately  8.5  km  MSL 
(Table  3).  The  wavelength  and  wave  amplitude  (i.e.,  in 
terms  of  w)  are  comparable  to  those  observed  in  IOP  4a. 
For  the  three  legs  between  4.5  and  6  km  MSL,  clusters  of 
sporadic  perturbations  are  evident  in  the  second  wave 
crest,  which  is  centered  over  Owens  Valley  and  has  sub¬ 
stantially  larger  amplitude  in  terms  of  the  vertical  dis¬ 
placement  than  the  first  and  third  wave  crests  (Fig.  14). 
Again,  this  suggests  that  wave-induced  KH  instability 
only  occurs  when  the  wave  amplitude  exceeds  a  certain 
threshold.  Similar  to  IOP  4a,  the  KH  at  the  highest 
altitude  leg  is  also  characterized  by  relatively  smooth 
kilometer-scale  eddies  separated  by  narrow  turbulent 
zones,  probably  resulting  from  the  intersection  of  the 
flight  segment  with  the  undular  top  of  breaking  KH  bil¬ 
lows.  Along  the  two  lower  KH  legs,  kilometer-scale 
eddies  are  evident  within  the  second  wave  crest,  em¬ 
bedded  in  strong  turbulence.  The  wavelet  spectra  of  the 
three  KH  legs  for  IOP  la  also  exhibit  the  vertically  co¬ 
herent  conelike  structures  with  distinctive  maxima  in 
period.  These  “cones”  point  toward  smaller  scales,  sug¬ 
gestive  of  a  predominant  downscale  energy  cascade  (Fig. 
15).  Small-scale  turbulence  and  large  eddies  are  evident  in 
the  wavelet  spectra  of  the  turbulence  leg  below  the 
mountaintop  and  these  large  eddies  are  primarily  shear- 
induced,  probably  associated  with  the  interaction  between 
waves  and  in-valley  flows.  It  is  evident  that  the  nonlinear 
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Fig.  13.  (a)  Profiles  of  the  cross-valley  and  along-valley  wind 
components  (i.e.,  U  and  V)  and  potential  temperature  derived  from 
the  1800  UTC  2  Mar  2006  NCAR  ISS  sounding  and  UWKA  mea¬ 
surements.  (b)  The  buoyancy  frequency  squared,  Scorer  parameter, 
and  gradient  Richardson  number  derived  from  the  sixth-order 
polynomial  regressions. 

index  Nrj/U  becomes  significantly  larger  below  the 
mountaintop  level  where  the  cross-valley  winds  are 
much  weaker  (Table  3).  As  demonstrated  in  previous 
studies,  for  highly  nonlinear  waves,  turbulence  could 
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Fig.  14.  The  25-Hz  vertical  velocity  for  IOP  la  is  plotted  as 
a  function  of  the  horizontal  distance  from  a  reference  point  over 
the  Sierra  ridge  (the  west  end-point  of  the  top  leg)  for  legs  (top)- 
(bottom)  1-9  and  17,  including  topography.  The  corresponding 
vertical  displacement  and  filtered  vertical  velocity  curves  are 
shown  as  dashed  and  gray  curves,  respectively. 


extract  energy  from  waves  associated  with  nonlinear 
processes  such  as  wave  steepening  or  breaking  (Fua  et  al. 
1982;  Weinstock  1986).  The  interaction  between  trapped 
waves  and  a  stagnant  layer  beneath  a  stable  shear  layer 
has  been  examined  by  Jiang  et  al.  (2006).  Using  a  me- 
soscale  model,  they  demonstrated  that  turbulence 
develops  in  an  initially  nonturbulent  shear  or  stagnant 
layer  beneath  the  first  few  wave  crests,  and  the  trapped 
waves  decay  with  downstream  distance. 

For  the  IOP  la  flight,  the  third-order  structure  func¬ 
tions  shown  in  Fig.  16  suggest  a  possible  upscale  energy 
transfer  along  both  the  KH  and  in-valley  turbulence  legs 
with  horizontal  scales  ranging  from  a  few  hundred  me¬ 
ters  to  a  few  kilometers.  The  horizontal  scale  ranges  in 
which  the  upscale  energy  transfer  may  occur  are  in¬ 
cluded  in  Table  3.  Negative  dissipation  rates  or  upscale 
energy  transfer  have  been  documented  by  previous  studies 
for  scales  between  500  and  1000  km,  which  has  been  in¬ 
terpreted  using  the  two-dimensional  turbulence  theory 
or  attributed  to  the  Coriolis  effect  (Hogstrom  et  al.  1999; 
Cho  et  al.  2001).  More  recently,  Koch  et  al.  (2005)  an¬ 
alyzed  aircraft  measurements  of  turbulence  associated 
with  gravity  waves  excited  by  a  frontal  system.  Negative 
energy  dissipation  was  found  for  horizontal  scales  from 


Table  3.  As  in  Table  1,  but  for  IOP  la. 


Leg 

Time  (UTC) 

Alt  (m) 

Char 

Data 

points 

a 

Upscale 
range  (m) 

1 

1659-1708 

8550 

Wave 

14  275 

0.25 

450  up 

2 

1711-1724 

7640 

Wave 

19  025 

0.23 

430-830 

3 

1726-1734 

7050 

Wave 

12  025 

0.23 

— 

4 

1736-1750 

6110 

Wave 

21  025 

0.23 

— 

5 

1753-1800 

5200 

KH 

11  275 

0.27 

530  up 

6 

1803-1814 

4870 

KH 

17  275 

0.26 

— 

7 

1815-1820 

4570 

KH 

8275 

0.29 

— 

8 

1824-1834 

4280 

Turb 

15  525 

0.58 

— 

9 

1835-1840 

3960 

Turb 

7525 

1.3 

— 

10 

1841-1847 

3660 

Turb 

9775 

2.5 

— 

11 

1848-1851 

3350 

Turb 

4525 

— 

300-730 

12 

1852-1856 

3040 

Turb 

6025 

— 

3500  up 

13 

1857-1900 

2710 

Turb 

3775 

— 

200  up 

14 

1901-1904 

2460 

Turb 

5275 

— 

— 

15 

1905-1908 

2140 

Turb 

4775 

— 

1700-2200 

16 

1909-1912 

1840 

Turb 

4525 

— 

250  up 

17 

1913-1916 

1450 

Turb 

4525 

— 

280  up 

300  to  700  m,  although  no  physical  explanation  was  pro¬ 
posed.  According  to  the  wavelet  analysis,  the  downscale 
energy  cascade  associated  with  KH  instability  exhibits 
scale  preferences  and  generates  distinctive  energy  max¬ 
ima  at  different  scales.  The  energy  cascade  associated  with 
steepening  or  breaking  mountain  waves  may  have  similar 
characteristics  (Weinstock  1986;  Jiang  and  Doyle  2004; 
Doyle  et  al.  2005).  In  contrast  to  classical  phenomeno¬ 
logical  turbulence  theory,  which  assumes  steady  forcing 
(i.e.,  energy  source)  in  the  large-  or  small-scale  limit,  the 
energy  maxima  associated  with  KH-  or  wave  breaking¬ 
generated  large  eddies  may  serve  as  transient  energy 
sources  and  disperse  energy  toward  both  larger  and  smaller 
scales. 

7.  Summary 

The  in  situ  measurements  from  two  UWKA  flights  over 
Owens  Valley  during  the  Terrain-Induced  Rotor  Exper¬ 
iment  have  been  analyzed  with  a  focus  on  wave-induced 
KH  instability  and  turbulence  characteristics  in  the  lower 
to  middle  troposphere  over  complex  terrain.  The  two 
events  examined  are  characterized  by  interaction  be¬ 
tween  a  stably  stratified  shear  layer  above  the  moun- 
taintop  level  and  stationary  moderate-amplitude  trapped 
waves.  According  to  the  observations  and  linear  analysis, 
KH  instability  develops  in  the  crests  of  trapped  waves 
where  the  Richardson  number  is  substantially  reduced. 
KH-like  eddies  are  absent  in  wave  troughs,  where  the 
local  Richardson  number  is  increased.  Wave-induced  KH 
instability  also  shows  wave  amplitude  dependence.  In 
both  IOPs,  the  KH  eddies  are  observed  in  the  wave 
crest  with  the  largest  wave  amplitude,  implying  that  KH 
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Fig.  15.  As  in  Fig.  8,  but  for  the  three  KH  legs  and  a  turbulence 
leg  (i.e.,  leg  10)  of  IOP  la.  The  locations  of  the  second  and  third 
wave  crests  are  highlighted  by  dashed  ellipses. 

instability  tends  to  develop  only  when  the  wave  ampli¬ 
tude  exceeds  a  certain  threshold. 

The  phase  relations  between  different  velocity  com¬ 
ponents  and  potential  temperature  for  the  observed 
kilometer-scale  eddies  are  consistent  with  growing  KH 
billows,  suggesting  that  these  eddies  could  extract  energy 
directly  from  the  sheared  mean  flow.  The  coexistence  of 
large  eddies  and  strong  turbulence  observed  along  some 
of  the  KH  legs  suggests  those  KH  billows  are  likely 
breaking.  The  energy  cascade  associated  with  the  break¬ 
ing  KH  billows  exhibits  scale  preferences  and  results 
in  distinctive  maxima  in  period  (or  wavelength)  in  the 
Fourier  power  spectra  and  wavelet  spectra.  For  KH  legs 
characterized  by  small-scale  turbulence  embedded  in 
kilometer-scale  eddies,  the  mean  spectral  slope  becomes 
— 5/3  for  scales  up  to  1-2  km,  which  is  consistent  with  the 
turbulence  inertial  subrange  scaling.  It  is  noteworthy 
that  the  corresponding  second-  and  third-order  absolute 
structure  functions  are  also  consistent  with  the  inertial 
subrange  scaling  over  a  similar  separation  distance.  The 


r  (m) 

Fig.  16.  As  in  Fig.  11,  but  for  legs  (a)  5  and  (b)  16  of  IOP  la. 

analysis  of  third-order  structure  functions  suggests  a 
possible  upscale  energy  transfer  on  scales  from  a  few 
hundred  meters  to  a  few  kilometers,  possibly  associated 
with  the  transient  energy  dispersion  from  breaking  large 
eddies  generated  by  KH  instability.  These  KH  eddies 
are  generally  anisotropic  with  the  vertical  to  horizontal 
aspect  ratio  less  than  unity  because  of  the  stratification 
effect.  They  are  characterized  by  large  negative  skewness, 
indicative  of  strong  asymmetry  in  updrafts  and  down- 
drafts  associated  with  the  positive  vertical  wind  shear. 
It  is  noteworthy  that  large  eddies  were  also  observed 
below  the  mountaintop  level  during  T-REX  IOP  13 
(Doyle  et  al.  2009).  Those  eddies,  or  subrotors,  were 
induced  by  shear  instability  associated  with  an  elevated 
horizontal  vortex  sheet  generated  through  boundary  layer 
separation  (Doyle  et  al.  2009). 

Below  the  mountain  top  level,  vigorous  turbulence  is 
probably  developing  because  of  interaction  with  trap¬ 
ped  waves  aloft,  which  may  have  caused  the  fast  decay 
of  the  trapped  waves  downstream.  Turbulence  extracts 
energy  from  waves  through  nonlinear  wave  steepening 
or  wave  breaking  processes  in  accordance  with  the  rel¬ 
atively  weak  cross-valley  winds  below  the  mountain  top 
level.  The  associated  energy  cascade  is  predominantly 
downscale,  but  with  scale  preferences  that  result  in  a 
number  of  kilometer-scale  energy-containing  maxima  in 
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the  power  spectra  and  wavelet  spectra.  Similar  to  KH 
eddies,  these  distinctive  energy  maxima  may  serve  as 
transient  energy  sources  and  transfer  energy  toward  both 
larger  and  smaller  scales.  In  general,  the  turbulence  be¬ 
low  the  mountaintop  level  is  better  developed  and  more 
isotropic  than  that  at  the  KH  levels,  and  it  generally  fol¬ 
lows  the  —  5/3  power  law  for  scales  less  than  1  km. 
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Fig.  Al.  The  local  Richardson  number  computed  using  (A4)  is 
contoured.  The  contour  interval  is  0.1.  Relevant  parameters  used  in 
this  example  can  be  found  in  the  appendix. 


APPENDIX 


Ri(x,  z)  =  Ri. 


0,  G  8 


(A4) 


Richardson  Number  Modified  by  a  Trapped  Wave 

Consider  the  interaction  between  a  shear  layer  char¬ 
acterized  by  a  constant  buoyancy  frequency  N  and  a  lin¬ 
ear  shear  Uz  =  constant  and  a  simple  stationary  wave 
trapped  between  two  rigid  lids  located  at  mz  =  ±tt!2  in  an 
x-z  plane.  The  vertical  velocity  induced  by  the  trapped 
wave  can  be  written  as  w'(x,  z)  =  w0  cos  (mz)  cos  (kx), 
where  m  and  k  are  the  vertical  and  horizontal  wave- 
numbers,  respectively,  and  wq  is  the  maximum  vertical 
motion.  The  ambient  Richardson  number  is  Ri0  =  N2/U2. 
Using  linear  theory,  we  obtain  the  horizontal  velocity  and 
potential  temperature  perturbations  u'  and  0'  and  the 
vertical  displacement  77: 

r)(x,  z)  =  ^  cos  (mz)  sin  (kx),  (Al) 

mw 

u'(x ,  z)  =  — - —  sin(raz)  sin(kx),  and  (A2) 

K 

0f(x,  z)  =  — jjf  -  cos  (mz)  sin(kx).  (A3) 

U  rC 

Assuming  0'/0z,  2 u’JU v  g'/g  <  1,  where  g'  =  Uwx  de¬ 
notes  the  wave-induced  vertical  acceleration,  to  the  first- 
order  approximation,  the  local  gradient  Richardson  number 
can  be  written  as 


The  last  three  terms  in  the  brackets  correspond  to  wave- 
induced  stability  variation,  shear  variation,  and  wave- 
induced  vertical  acceleration  effects,  respectively.  Using 
(A1)-(A4),  we  obtain 


Ri(x,  z)  =  Rio+Rioi? 

2m2  U 


-jj  +  m  tan  (mz) 

&2c/2~ 

g 


(A5) 


Equation  (A5)  indicates  that  both  the  wave-induced 
shear  and  vertical  acceleration  (i.e.,  last  two  terms  in  the 
brackets)  tend  to  reduce  the  local  Richardson  number  in 
a  wave  crest  (i.e.,  77  >  0)  and  act  in  the  opposite  sense  in 
a  wave  trough.  In  addition,  the  Richardson  number 
variation  is  proportional  to  the  wave  amplitude. 

As  an  example,  we  use  a  set  of  parameters  estimated 
from  IOP  4a:  w0  =  6  m  s-1,  Uz  =  0.013  s-1,  k  =  3.14  X 
10“4  m"1,  U(0)  =  20  m  s"1,  and  N  =  0.013  s“\  which 
yields  Ri0  =  1.  The  local  Richardson  number  computed 
using  (A5)  and  m  =  yj ( N/U )2  —  k 2  is  shown  in  Fig.  Al. 
It  is  evident  that  the  Richardson  number  is  substantially 
reduced  within  the  wave  crest  and  the  lower  half  of  the 
wave  (i.e.,  z  <  0),  and  Ri  is  increased  in  wave  troughs.  The 
Richardson  number  minimum  is  located  at  77  =  77max  and 
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mz  ~  — 7t/4.  The  relative  importance  of  wave-induced 
shear  to  acceleration  can  be  expressed  as  m2g/(k2UUz). 
For  the  given  example,  the  acceleration  term  is  much 
smaller  than  the  other  three  terms. 
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